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Abstract 

We study the matching of LRS spatially homogeneous collapsing dust space-times with 
non-stationary vacuum exteriors in cylindrical symmetry. Given an interior with diagonal 
metric we prove existence and uniqueness results for the exterior. The matched solutions 
contain trapped surfaces, singularities and Cauchy horizons. The solutions cannot be asymp- 
totically flat and we present evidence that they are singular on the Cauchy horizons. 
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1 Introduction 

In cylindrical symmetry, unlike spherical symmetry, there exist exact solutions of the Einstein 
field equations (EFEs) which contain gravitational radiation. This raises the possibility of con- 
structing collapsing matter cylinders matched to vacuum but radiating exteriors. In particular 
one can attempt to find analogues in cylindrical symmetry of the Oppenheimer-Snyder model. 
The Oppenheimer-Snyder model results from the matching of a ball in the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) dust space-time to a Schwarzschild exterior across a timelike surface 
preserving the spherical symmetry. Thus there is no radiation in the exterior. This model 
cannot be generalised (in a physically reasonable way) to cylindrical symmetry with spatially 
homogeneous interiors and static exteriors. In fact, Mena, Tavakol and Vera [8] have shown that 
the matching of such space-times leads to very restrictive interior solutions with no evolution 
along one spacelike direction. As a next step one therefore considers a non-static exterior, but 
for simplicity retains a spatially homogeneous interior. 

If such a model is to represent a cylinder of dust collapsing in a vacuum, then it should 
be asymptotically flat. Berger, Chrusciel and Moncrief [1] have shown that for cylindrically 
symmetric vacuum spacetimes the requirement of asymptotic flatness implies that the 2-surfaces 
of transitivity are never trapped (specifically they show that one can choose a comoving radius 
of the form R(T, p) = p, which forbids trapping). 

In this paper we investigate a family of models which result from the matching of a spatially 
homogeneous collapsing dust interior with a non-stationary vacuum exterior. These interiors 
always contain cylinders of symmetry (which we shall call just 'cylinders' for brevity) which 
are trapped near to the singularity. The marginally-trapped cylinders trace out a 3-surface in 
the interior which eventually arrives at the boundary of the matter. Matching conditions at 
the boundary mean that the intersection 2-surface at the boundary is marginally-trapped as 



seen from the exterior and the results of Berger et al. [1] then prevent the exterior from being 
asymptotically flat. We shall argue that the failure of asymptotic flatness is in fact due to 
singularities in the vacuum exterior. 

The technique will be to use the interior, which can be assumed to be known, to give data 
on the boundary of the matter, which will be a co-moving cylinder from the point of view of the 
interior but just a time-like surface made up of cylinders of transitivity seen from the outside. 
The EFEs in the exterior are a system of hyperbolic equations in the 2-dimensional quotient 
space Q of the space-time by the symmetries (this is the (T, p)-space) and data can be given 
on any non-characteristic curve in Q. The boundary of the matter defines such a curve, even 
though it is a time-like surface in the full 4-metric. Therefore we can deduce existence and 
uniqueness of solutions in the domain of dependence in Q of the curve on which data is given. 
If the interior collapses to a singularity, as it does in these examples, then the data-curve has 
an end and there is a Cauchy horizon in the exterior. One expects the exterior to be singular 
on this Cauchy horizon, and we present some evidence for this belief. We shall interpret the 
singularity at the Cauchy horizon in terms of incoming gravitational radiation. 

The plan of the paper is as follows. In Section 2 we review the matching procedure in 
general, and in Section 3 we introduce the interior and exterior metrics which we wish to match. 
The matching conditions for these metrics are written out in Section 4. In Section 5, we give 
some properties of the matched space-times when the metrics are diagonal in the canonical 
coordinates, obtaining existence from standard theorems for linear PDEs. We briefly discuss a 
non-diagonal example in Section 6 and finally in Section 7 we discuss some properties of the 
solutions, arguing that the Cauchy horizon is singular and giving a physical reason for this. 
Some equations are collected into an Appendix. 

We use units such that 8irG = c = 1. Greek indices take values 1 to 3 and latin indices to 

3. 

2 Matching procedure in brief 

In this section, we recall how two space-times can be matched across a hypersurface by showing 
which junction conditions the respective metrics should satisfy. The first set of junction condi- 
tions ensures the continuity of the metrics across the matching hypersurface. The second set of 
junction conditions ensures the continuity of the second fundamental forms and, given the first 
set, is equivalent to requiring a non-singular Riemann tensor distribution, which prevents infinite 
discontinuities of matter and curvature across the matching hypersurface. In what follows we 
assume that the matching surface is non-null. For cases where the matching surfaces can change 
character see [7]. 

Let (V + ,g + ) and (V~,g~) be two C 3 space-times with oriented boundaries a + and a~ , 
respectively, such that a + and a~ are diffeomorphic. The matched space-time (V,g) is the 
disjoint union of V^ 1 with the points in identified such that the junction conditions are 
satisfied (Israel [4], Clarke & Dray [3] and Mars Sz Senovilla [7]). Since are diffeomorphic, 
one can then view those boundaries as diffeomorphic to a 3-dimensional oriented manifold a 
which can be embedded in V + and V - . Let and {x ±% } be coordinate systems on a and V ± 
respectively. The two embeddings are given by the following C 3 maps 

: a — ► V ± (1) 

r -> x i± = ^ ± (n, 

such that = $ ± (cr) C V^. The diffeomorphism from a + to o~ is <I> o 1 . 



Given the basis {d/d^ a \ p } of the tangent plane T p a at some p € a, the push- forwards d^\ p 
map {d/d£ a \ p } into three linearly independent vectors at $ ± (p) represented by e^\q,±( p y. 
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On the other hand, using the pull-backs of the maps the metrics can be mapped to 
a given two symmetric 2-covariant tensors g + and g~ whose components in the basis {d£ a } are 



9 a/3 — e t le p j 9ij\ a ± 
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The first matching conditions are given by the equality of the first fundamental forms (Israel 
[4]) " 

dtp = 9 af s- (2) 

We note that the existence of a continuous metric allows for the treatment of Einstein's field 
equations in the distributional sense (see e.g. Mars & Senovilla [7]). The bases {e^ \ p } and 
{e^|p} can be identified, 

= d*~ f 3 
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as can the hypersurfaces a + = a , so henceforth we represent both by a. 
We now define a 1-form n, normal to the hypersurface a, as 

n^e ±) = 0. 

The vectors {n ± ,e^} constitute a basis on the tangent spaces to at a^. Since the first 
junction conditions allows the identification of {e^} with {e^}, we only have to ensure that 
both bases have the same orientation and that nfn +t = n^n~ l is satisfied in order to identify 
the whole 4-dimensional tangent spaces of at a, {n + , e+} = {ft ~, e~}. 
The second fundamental forms are given by 



H 



aj3 



-n ± p ± i\7 ± p ±i 



and the second matching conditions, for non-null surfaces, are the equality of the second funda- 
mental forms 

Kp = H-p- (3) 
We note that these matching conditions do not depend on the choice of the normal vectors. 



3 Cylindrically symmetric exterior and spatially homogeneous 
interior 

We are interested in the particular cases for which the matching surface inherits a certain 
symmetry of the two space-times (V^, g^). Such matching is said to preserve the symmetry. In 
practice one demands that the matching hypersurface is tangent to the orbits of the symmetry 
group to be preserved (see e.g. [11]). In the case we shall consider below, the symmetry group 
is an Abelian G2 and we are interested in the matching of a cylindrical symmetric exterior to a 
spatially homogeneous (and anisotropic) interior preserving the G2 symmetry. 

Since we want to preserve the symmetry of a G2 exterior space-time we must ensure that 
the interior also has the same group of symmetries, possibly as a subgroup of a larger isometry 
group. We wish to preserve an axial symmetry so that we require the interior to have an axial 



symmetry in the sense of a periodic Killing vector, vanishing on a symmetry axis. In the exterior 
V + the Killing vectors commute, so we require this in the interior. Thus we require the interior 
to be spatially homogeneous, with an axial symmetry and at least one Killing vector commuting 
with the axial symmetry. This is sufficient to force the interior to be an LRS space-time (see 
Mena, Tavakol & Vera [8]). 

We deal separately the metric forms for the interior and exterior. 



3.1 Interior 

All the LRS spatially homogeneous metrics can be written in the compact form 
ds 2 ~ = -dt 2 + a{t) 2 2 + b{t) 2 [{dr - erdzf + E(r) 2 d^ 

where 

6 = dz + n(F(r) + k)d<p 
and the functions E and F are given by 



(4) 



E(r) 



sinr, k = +1 

r, k = 
sinhr, k = — 1 



and F(r) 



— cos r, k = +1 
r 2 /2, k = 
coshr, k = — 1, 



and where e and n are given such that 

e = 0, 1; n = 0, 1; en = ek = 0. 

We note that 



(5) 



E = F r ; (E r ) 2 + /cE 2 = 1, 

where a subscript denotes the partial derivative with respect to the indicated variable. The axial 
Killing vector is then given by 

m = 9<p, 

while the other three Killing vectors ffi, i = 2 to 4 are 

m = d z , 

ffy = smipe tz d r + cosip(e ez f(r)d v + g(r)d z ), 
f]4 = cos(pe tz d r - sin (p(e ez f(r)d v + g(r)d z ), 

where we have defined f(r) = E r /E and g(r) = n(E — f(F + k)). The Killing vector which 
commutes with the axial Killing vector 771 is 772- We now assume a perfect fluid interior and 
recall a result which is a direct consequence of the so-called Israel conditions n~ l T^ = n +l T^: 



Lemma 1 Let (V,g) be a space-time resulting from the matching of two space-times (V ± ,g ,± ). If 
(V~, g~) has a perfect fluid and (V + ,g + ) vacuum then the perfect fluid has p = 0. Furthermore, 
if(V~,g~) is spatially homogeneous then p = everywhere. 

The metric for the interior will be assumed to be given by (4) with dust as matter content. For 
convenience we shall use the form 



ds 



2- 



-A 2 dt 2 + B 2 dr 2 - 2erB 2 drdz + C 2 dip 2 + 2Edzd<p + D 2 dz 2 , 



(6) 



KS/Bianchi types 


e 


n 


k 


I 











KS 








1 


KS,III 








-1 


IX 





1 


1 


II 





1 





VIII,III 





1 


-1 


v,vn fc 


1 









Table 1: Classification of the possible G3 on S3 subgroups according to the values of {e,k,n} 
for the metric given by (4). KS denotes Kantowski-Sachs. 



where A, B,C, D and E functions of t and r. The line-element (4) is recovered by making the 
identifications 



A 2 (t,r 
B 2 (t,r 
C 2 (t,r 
D 2 (t,r 
E(t,r 



1, 

b 2 (t), 

b 2 (t)Z 2 (r, k) + na 2 (t)(F(r, k) + k) 2 , 
a 2 {t) + er 2 b 2 {t), 
na 2 (t)(F(r, k) + k). 



In the following we shall take the functions in (6) to be arbitrary functions, with 

eE = 0, 

which follows from (5). 
3.2 Exterior 

The exterior will be assumed to be cylindrically symmetric vacuum (with an Abelian G2 on S2). 
The most general metric form can be written as [1] 



ds 2+ = e 2 ^\-dT 2 + dp 2 ) + R 2 e-Wd<p 2 + e 2 ^(dz + Wdtpf 



(7) 



where ifi, 7, R, W are functions of the coordinates p, T. For W 7^ the two Killing vectors are not 
hypersurface orthogonal and the cylindrical gravitational waves have two polarisations states. 
One of the EFEs is (33) 

Rtt — Rpp = 0, 



which has the general solution 
where 



R(T,p) = F(x) + G(y), 



(8) 



x = T + p; y = T - p. 
The metric (7) is invariant under a coordinate transformation 

t : (x,y) -» (f(x),g(y)), 



(9) 



where / and g are arbitrary differentiate functions (with non-zero derivative). One may use 
this freedom to restrict R, but we shall make a different choice below, which is to use r to 



prescribe the matching surface a + . The choices W = and R = p in (7) lead to the well-known 
Einstein-Rosen metric. 

It will be important below to see when the 2-surfaces of transitivity are trapped or marginally- 
trapped. From (7) we see that these 2-surfaces have intrinsic metric Rdipdz. Taking the two 
null-normals to these 2-surfaces, we see that such a 2-surface is trapped if 

AR x Ry = (R T - R P )(R T + R P ) > 0, (10) 

and is marginally-trapped if this expression is exactly zero. It follows from (8) that if R x = 0, 
at say x = xo, then all 2-surfaces of transitivity with x = xq are marginally-trapped. Similar 
statements follow with y replacing x. 



3.3 The matching hypersurfaces 

The embeddings a ± can be defined by using coordinates on a denoted by £ a = {\,4>, (}■ We 
choose (j) such that 











^ \d<p) dp 

Since we perform the matching preserving the Gi symmetry there must a vector field 7 which 
together with d/dp\ a - generates the G2 on S^- In particular 7 must satisfy 



d& (7) = ci 
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with ci,C2 constants such that C2 / and otherwise arbitrary. We now use the fact that d/d(j) 
and 7 must commute to choose 7 = d/d( and a coordinate transformation (' = C2C, 4>' = 4> + ciC 
to get (after dropping the primes): 
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Finally, A is chosen such that d& (d/d\) is orthogonal to e 2 and e 3 . This implies 
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By denoting the embedding {$° , <I> 1 , <I> 2 , <I> 3 } in (1) as {t, r, ip, z}, the matching surface a 
is parametrized as 

= r ^> z) ■ t = t{\), r = r(A), p = 4>,z = C + fz(X)}, 

where t(X) and r(A) are functions of A restricted by the fact that d<f>~ has to be of rank 3, that 
is i 2 + r 2 ^ 0, and 

A(A) ^ er(A)r(A)^ (11) 

where the dot denotes differentiation with respect to A. 

Consider now the embedding <E> + . The axial killing vectors from both must coincide at 
the matching surface so 











and d& + (d/d()\ a+ must complete the basis of the G2 on S2. By a similar procedure as for a 
we can always choose C and A such that 



d_ 

d(f> 



d_ 



+ dX dp~ 



<T+ 



dX dT 

so that matching surface a + can be parametrized by 

a+ = {( T > P, <M) : T = T W,P = ^ = </>, z = C}. 



4 The matching conditions 

In order to derive the junction conditions we have to calculate the first and second funda- 
mental forms for both a + and o~ . For the g~ metric, the parametric form of o~ gives 
dt\ a - = idX, dr\ a - = rdX, dip\ a - = dcf) and dz\ a - = d( + f z dX. Using (11), the first funda- 
mental form g~ b on u~ can be written as 

ds ^-\ a _ 1 (_^ 2 i 2 + B 2 r 2 )dX 2 + C 2 d<f + 2Ed<t>d(; + D 2 d( 2 , 

where 

B 2 ^B 2 (l-^), 

Similarly the first fundamental form g + on a + is given by 

ds 2+ \ a+ £ e 2 ^\-f 2 + p 2 )dX 2 + (R 2 e-^ + IU 2 e 2 ^)# 2 + 2We 2 ^d(d<P + e 2 ^dQ 2 

and then the first set of matching conditions (2) is: 

-A 2 i 2 + B 2 r 2 ^ e 2 ^\-f 2 + p 2 ) (12) 

C 2 £ R 2 e- 2 ^ + W 2 e 2 ^ (13) 

D 2 ± e 2 ^ (14) 

E = We 2 ^ (15) 

The normal forms to the matching surfaces can be taken as 

n = AB(-rdt + idr)\ a - 

n + = e ^-Tp)(-pdT + Tdp)\ a+ 

In our case, the interior space-time contains dust and so the matching is performed across a 
time-like hypersurface ruled by matter trajectories, which are geodesies (this follows from the 
matching conditions). Therefore, t = 1 and r = 0. From the point of view of the exterior, the 
matching surface is given parametrically as (x(X),y(X)) in terms of the coordinates x,y in (9). 
The product xy does not vanish since the matching surface is everywhere time-like and so we 
may exploit the transformation r of (9) to set x = A, y = X, or equivalently T = 1 and p = 0. 
Now the matching surface is at r = ro, p = po and the time-coordinates can be taken to agree, 
t = T. 



Under these circumstances we get the first set of matching conditions (12)-(15) in the form: 



7 = i> (16) 

b 2 ^ 2 + na 2 {F + k) 2 1 Ri e -W + w 2 e 2 <C (17) 

a 2 + er 2 6 2 a ^ (lg) 

na 2 (F + k) = We 2 ^, (19) 



and the second set of matching conditions (3) as 



% = (20) 
7p = (21) 



na 2 F, 



r 



b 



and for (DB) t + 



e 2 ^W p (22) 
RRpe- 2 * + WW p e 2 ^, (23) 



As a consequence of (16) and T = 1 we also have 

7t = IpT- 

The fact that e = excludes some Bianchi types: 

Proposition 1 A G4 on S3 LRS spacetime admitting a simply transitive subgroup G3 of Bianchi 
types V or Vllh cannot be matched to a cylindrically symmetric non- stationary vacuum space- 
time across timelike surfaces preserving the cylindrical symmetry. 

We shall now look more closely at the cases with diagonal metrics. 



5 Diagonal cases 

For simplicity, we begin with a diagonal interior metric, so that e = n = 0. Now W and W p 
vanish on a, from (19) and (22). In each case below we first establish that R 7^ in the exterior, 
and we can then deduce that W = in the exterior, from uniqueness for the linear equation 
(35) which it satisfies. Assuming this for now, the matching conditions give 

e^ = a (24) 
R = afeS (25) 
R p = aS r , (26) 

together with 7 = ip and 7^ = ip p = 0. We first consider the two explicit solutions corresponding 
to the flat (k = 0) FLRW and to the Bianchi I cases and then consider implicit solutions for 
Kantowski-Sachs and Bianchi III models. We may summarise the results that we shall find in 
the following theorem. 

Theorem 1 A cylindrical interior dust metric constructed from a collapsing FLRW, Bianchi 
I or Kantowski-Sachs cosmological model can be matched to a diagonal vacuum exterior within 
the domain of dependence V of the matching surface. The metric is smooth up to the Cauchy 
horizon TL, but there are trapped surfaces in the exterior. If the interior is a Bianchi III model 



(with space sections IR x S 2 ) then the matching can be performed as far as the Cauchy horizons 
of the past and future singularities provided the matter occupies less than half of the S 2 . If 
the matter occupies more than half of the S 2 then the matching is at least possible near to the 
matching surface. 

In Section 7 we shall present evidence that the exterior is in fact singular on the Cauchy surface, 
and we shall suggest a physical origin for this singularity. 

5.1 FLRW k = 

We want the interior metric to be collapsing so that the relevant solution to EFEs (38) is 
explicitly given by 

a(t) = b(t) = (a - i) 2/3 , (27) 

for t G (— oo,a). The matching surface a is the cylinder p = po, T < a, terminating in the 
singularity at T = a. At a, by (25) and (26) we have 

R g= ro (Q,_T) 4 / 3 (28) 
R p ^ (q-T)I, (29) 

so that in particular R and R p on a are positive for T < a, vanishing only at T = a. 

We shall determine the exterior using the matching conditions, in the domain of dependence 
T> = {T + p < a + po, p > po} of a in Q, when a is taken to be at p = po- The boundary of 
the domain of dependence, at T + p = a + po, is the Cauchy horizon which we shall call Ji. It 
follows from the D'Alembert solution of the one-dimensional wave equation for R that, since the 
data for it is positive, R is positive in V. (Now we may conclude that W = 0). We may obtain 
R explicitly from (28), (29) and (8) as 

R=^(a + po-T-p)l-^(a + po-T-p)l + r ^(a + po-T + p)l+^(a + po-T + p)l. (30) 

From this formula for R we can investigate the condition (10) for trapped-ness. On a it is easy to 

/„ \3/2 

see that R y vanishes only at T = a, while R x vanishes at T = a and at T := To = a — (^f 1 ) 
Thus there is a marginally-trapped cylinder at T = To on a and this gives rise in the exterior 
to whole null hypersurface Ho of marginally-trapped cylinders at T + p = Tq + po- Cylinders to 
the past of Mo are not trapped and those to the future are. 

We have seen that R is non-negative in V and by inspection it vanishes only in the corner 
at T = a, p = po, and is smooth in the interior of T>, though not along the Cauchy horizon TL at 
T + p = a + po where second and higher derivatives diverge. 

At a, by (20), (24) and (27) we have 

V> = fln(a-T) 
^ = 0. 

This is the data for ■0 in the exterior, and it is smooth (in fact analytic) up to T = a. Given 
what we know about R, it now follows from standard theorems for linear hyperbolic equations 
(see e.g. section 6.5 in [10]) that a unique Y> exists satisying (34) in the interior of V, with 
the given data on a. There is no reason to expect the solution to be well-defined at TL, since 
it diverges on a at T = a, and in fact we shall argue in section 7 that there is a curvature 
singularity along H. 



Since 7 = tp and 7 p = ip p , once we have tfj in V the same argument applied to (32) gives 
existence and uniqueness for 7 in V, and we have constructed a vacuum exterior which matches 
on to the dust interior, at least inside T>. (In this and all the examples, the constraint equations 
(36) and (37) are identically satisfied by virtue of the field equations in the interior, as must be 
the case.) 

5.2 Bianchi I 

The collapsing solution to the EFEs (38) in the interior can be taken to be 

a(t) = (a-t)- 1/3 (/3-t) 
b(t) = (a-t) 2/3 

for t < min{a,/3} and r < tq. If a = (3 this reduces to the previous case. The character 
of the solution depends on which is the smaller of a and (3. If a < (3 then a diverges and b 
vanishes at the singularity. This is the behaviour of the Kasner (vacuum) solution and we can 
call it Kasner-like. If a = (3 then a and b vanish at the same rate, which is precisely the FLRW 
behaviour so that we may call it FLRW-like. If a > (3 then a vanishes but b does not and we 
call this disk-like. 

Since the interior is given explicitly we can use the matching conditions as before to get the 
exterior functions. We again find that the data R and R p are positive so that R is positive in T> 
and W = 0. Explicitly, (25), (26), (8) imply 

R = ^(a + po-T - p)3(/? + po-T - p) -\{a + p - T - P )*(J3 + Po -T - p) 
+^(« + Po-T-/>)3 

+y (a - Po - T + p)\tf - po - T + p) + ^(a - p - T + p)%((3 - p Q - T + p) 
-2q(« + Po -r + p)s. 

It is straightforward to see that R y is always negative on a while R x necessarily changes sign. 
Thus, by (10), there is always a marginally-trapped cylinder on a and this gives rise to a null hy- 
persurface of marginally-trapped cylinders A/o as in the previous case (for suitable combinations 
of the parameters a, {3 and po there may be three such hypersurfaces). However R is positive in 
the exterior, vanishing only in the corner p = po, T = min(a,/3). Note also that R x diverges at 
the singularity, and therefore on the Cauchy horizon, if the singularity is Kasner-like. 
For the other exterior metric functions we get from (24) 

Y> = 7 = ln(/?-T)-^ln(a-T) 

and with 7 p = ip p = we again have existence and uniqueness of solutions in the interior of T>. 
As before, we expect the solution to diverge near Ti, a point we shall return to below. We have 
then constructed a vacuum exterior which matches on to the Bianchi I dust interior, at least 
inside V. 

5.3 Kantowski-Sachs and Bianchi III 

For these interiors, implicit solutions are known and we give them in the appendix. Dealing 
first with the Kantowski-Sachs metric, the data for R are positive on a so that as usual R is 



positive in the exterior, out to the Cauchy horizon, and W = 0. The singularity is Kasner-like, 
FLRW-like or disk-like according as Co is positive, zero or negative in (40). In each case, the 
quantity Rt — R p changes sign from positive in the infinite past to negative near the singularity 
so that there are trapped cylinders on a near the singularity and these persist in the exterior as 
before. The existence argument carries through just as before, out to the Cauchy horizon. 

There are two extra complications with the Bianchi III metric. The first is that there are 
singularities in the past as well as in the future. The second is that £ = sin r so that R p = a cos vq 
which can be negative or zero. 

Suppose the past and future singularities lie on a at T = Tq and T = T\ respectively, then 
the domain of dependence is restricted to 

V = {T + p < Ti + po, T - p > T - p } 

and there is a past Cauchy horizon as well as a future one. The nature of the singularities 
depends on the relative sizes of the constants M and Co in (39) and various combinations of 
Kasner-like, disk-like and FLRW-like are possible. However, at least one of the past and future 
singularities is always Kasner-like. 

Now for R, suppose first that < ro < 7r/2, then the data for R are positive so R is positive in 
the exterior and W = 0. The proof of existence is as before. On the other hand, if ir/2 < vq < ir 
then R p is negative on a and so R may vanish inside T>. If R vanished in a smooth way, this 
could just be the sign of a (second) regular axis for the axial Killing vector. To see that this 
is not the case we recall from the discussion of the Bianchi I case above that the derivative of 
R is singular on the Cauchy horizon emanating from a Kasner-like singularity. Since one of 
the singularities must be Kasner-like, R cannot be smooth at any axis meeting that part of the 
Cauchy horizon. 

For each kind of singularity there are both future-marginally-trapped and past-marginally- 
trapped cylinders on a. This follows from an analysis of the implicit solution (39), but can be 
seen just by considering asymptotic forms near the different kinds of singularity as they have 
already occurred. The global picture is complicated by the presence of both past and future 
Cauchy horizons, but the existence of the exterior in the whole domain of dependence follows 
as before if ro < 7r/2 and follows near to a if ro > tt/2. 

6 A Non-diagonal Example 

This case is much harder and we have only partial results. In principle, one could repeat the 
work of the previous section with the non-diagonal interior metrics of Bianchi types II, VIII 
and IX matched to non-diagonal exterior vacuum metrics. This is hampered by two things: the 
non-existence, in the literature at least, of even implicit solutions for the interior metric; and 
the fact that the EFEs are now coupled and cannot be solved one by one. 

For the interior solution, the Bianchi type II has been taken furthest and we shall consider 
just some particular instances of this case. The constants in (4) are e = k = and n = 1 so that 
the matching conditions from (16)- (23) are 
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together, as always, with 7 = ip and 7 p = ip p . 

The field equations for the interior metric are given by (41). They have been solved in terms 
of fi subject to a single ODE by Maartens and Nel [6] as follows: 

a(t) = /xexp(— 2 J n(s)(a — s)ds) 
b(t) = fi~ 1 exp(J n(s)(a — s)ds) 

where [i(t) satisfies 

Utt ~ -(fit? + 8(q - t)n/H - \ii 2 - 6(a - t) V = 0. (31) 

In (31), a is an additive constant and the singularity will be at t = a. It is possible to reduce 
the order of the equation for fi, but it does not seem possible to find a general solution. Two 
particular solutions are given by fi(t) = k/(a — t) 2 with k = 5/4 or k = 4/3, and now a and b 
are just powers of (a — t). 

Proceeding as before, we find that R is strictly positive, and in fact analytic, in the domain 
of dependence. There are trapped surfaces on a between the singularity at T = a and a value To 
which depends on the constants of integration. Thus there are trapped surfaces in the exterior. 

To obtain the other functions appearing in the exterior metric we need to solve the system of 
equations (32)- (37). This is a much-studied system and a good deal is known. Local existence 
in a closely related setting was shown by Chrusciel [2]. To obtain a more global result, and 
in particular existence out to the Cauchy horizon, the methods of Berger et al. [1] would be 
expected to work, but we have not done this. 



7 Properties of the solutions 

We shall concentrate on the diagonal cases, where existence is known. We have found that a 
radiating exterior can be matched to a collapsing interior but there are always trapped surfaces 
in the exterior. From the work of Berger et al. [1] this means that the solutions cannot be 
asymptotically flat, so that we have the following corollary to Theorem 1: 

Corollary 1 Asymptotically flat vacuum cylindrically symmetric spacetimes cannot be matched 
to the boundary of solid cylinders in LRS spatially homogeneous perfect fluid spacetimes with 
diagonal metrics, preserving the cylindrical symmetry. 

We now want to present some evidence that these solutions are actually singular on the Cauchy 
horizon H. We begin by introducing the quantity u defined by 

u = R 1 / 2 ^, 

then (34) becomes 

RxRy 

Uxy ~ AR 2 U ' 

It follows from this equation that if R x = 0, R 7^ at some value x = xq then u x is constant 
on x = xq, and therefore so is R l / 2 ipx- In the FLRW case, precisely this happens on H so 
that R 1 / 2 ^ x is constant on 7i. However this quantity diverges in the corner where H meets the 
matching surface, and it is therefore singular on H. Once ip is not differentiable at Ti, there is 
no reason to expect 7 to even be finite there, in which case the metric is singular. 



One can alternatively look at curvature components, specifically at the Weyl spinor in the 
exterior. We introduce the null tetrad 



m 



n 



I = 2- 1 / 2 e^(d T + d p ) 
i = 2- 1 ' 2 ^-\d T - d p ) 
i = 2- 1 / 2 (^d z + \^d <t> ), 



and calculate the (nonzero) components of the Weyl spinor in this tetrad as 
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If we restrict these quantities to a for the case of an FLRW interior, then we find that \I>rj = 
^4 = but that ^2 diverges as the singularity is approached. This is in the corner where H 
meets a and then we expect this singularity to propagate along the Cauchy horizon. 

We also note the appearance of the term R xx in ^q, which blows up at TC for all our examples. 
We cannot explicitly evaluate the other terms in \£o, but this suggests that \I/o diverges at H. 
This component of the Weyl spinor, in this tetrad, is interpreted as incoming radiation, which 
therefore becomes singular at the Cauchy horizon. This interpretation leads us to a physical 
explanation of the exterior singularity: suppose at a particular instant we set up Cauchy data 
with a spatially homogeneous interior, and we want the interior to remain spatially homogeneous 
even as it collapses. This requires incoming radiation and these examples show that it requires 
more and more as the interior approaches the singularity, leading to a singularity in the exterior 
on the Cauchy horizon. 
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Figure 1: Schematic diagram of the spacetime structure. 
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Appendix 

The EFEs for the cylindrically symmetric vacuum metric (7) are 



= lTT - 7pp -^ + ^ T -^-,((W T ) 2 -(W p ) 2 ) 



= Rtt-R pp 

= ^ TT + ^ T -^ pp -^ p -^{(W T f-(W p f) 

-r Wt + 1i 



= Wtt-W pp --!-W t + -?-W p + AW t1t -AW p1p 



together with the two constraint equations 



7p = R 2 ^ R 2 ( RR p(^T + ^p) ~ 2RRtM P + RpRpp - RtRtp) 

It = r2 \ R2 (RRt (V't + tf) ~ 2RR P ifrip P + RtR pp - R P R Tp ). 
The EFEs for LRS spatially homogeneous diagonal dust metrics (4) are 



fi = 


a t b t b 2 t k 
ab b 2 b 2 
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Jtt b 2 k 
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btt atbt (Ht_ 
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The solution for k = 1 (Bianchi III) in parametric form is [9] 



a cos r/ = M(r] sin rj + cos rj) + Co sin r\ 
b = ( 

where M and C\ are positive constants. 



b = C\ cos 2 rj, dt = 2bdr] 



The solution for k = — 1 (Kantowski-Sachs) in parametric form is [9] 

asinhry = [M(i] cosh rj — sinhr/) + Co cosh 77] 
b = C\ cosh 2 rj, dt = 2bdr] 

with M and C\ positive constants. 

The EFEs for LRS spatially homogeneous non-diagonal dust metrics (4) of Bianchi II 

b\ a t b t a 2 

^ ~ ~b 2+ ~aT~W 

btt b\ 3a 2 

n — 2— -i - 

U " b + b 2 46 4 

n + + + q2 

b a ab 46 4 



from which the solution of Maartens and Nel [6] given in the text can be found. 
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